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Non-equilibrium critical relaxation of the 3D Heisenberg magnets 
with long-range correlated disorder 

Pavel V. Prudnikov*), Maria A. Medvedeva, 
Dept. of Theoretical Physics, Omsk State University, Omsk 644077, Russia 

Monte Carlo simulations of the short-time dynamic behavior are reported for three- 
dimensional Heisenberg model with long-range correlated disorder at criticality, in the case 
corresponding to linear defects. The static and dynamic critical exponents are determined 
for systems starting from an ordered initial state. The obtained values of the exponents are 
in a good agreement with results of the field-theoretic description of the critical behavior of 
this model in the two-loop approximation. 

§1. Introduction 

Critical properties of disordered systems with short-range (SR) and long-range 
(LR) correlated randomness have been studied extensWely^^^^^^^^^^^ 
One important question to address is whether the introduction of weak randomness 
changes the universality class of phase transition. According to the well-known Har- 
ris criterion^* disorder with SR correlations is relevant if 2 — duQ = ao > 0, where d 
is the spatial dimension, and uq and ao are the correlation-length and the specific- 
heat exponents of the pure system. This criterion is modified in the presence of LR 
correlations in the disorder. A special type of such a disorder has been considered 
by Weinrib and Halperin (WH).^' They showed that the disorder with power law 
correlation g{x) ~ for large separations x is relevant if 2 — ovq > for a < d, 
whereas the usual SR Harris criterion recovers for a > d. As a result, the exis- 
tence of LR correlations in the disorder gives significant effect and wider class of 
disordered systems, not only the three-dimensional diluted Ising model with point- 
like uncorrelated defects, can be characterized by a new universality class of critical 
behavior. 

The power law decay for the impurity-impurity pair correlation function g{x) 
allows a direct geometrical interpretation. So, for integer a it corresponds to the 
lines (at a = d — 1) or the planes (at a = d — 2) of impurities of random orienta- 
tion-f^ Moreover, non-integer a is sometimes treated in terms of impurities fractal 
dimension.!^ Therefore, the models with LR-correlated quenched defects have both 
theoretical interest due to the possibility of predicting new types of critical behavior 
in disordered systems and experimental interest due to the possibility of realizing 
LR-correlated defects in the ^ He in aerogels,^ polymers,'^ and disordered solids 
containing fractal-like defects^^' or dislocations near the sample surface.'^ 

The critical exponents were calculated in Ref.^ in the one-loop approximation 
using a double expansion ine = 4 — (i<Cl and (5 = 4 — a ^ 1. The correlation- length 
exponent was evaluated in this linear approximation as = 2/a and it was argued 
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that this scahng relation is exact and also holds in higher order approximation. In 
papei^^ a renormalization analysis of WH model was carried out directly for the 
3D systems in the next two-loop approximation with the values of a in the range 
2 < a < 3. The static and dynamic critical exponents were calculated with the use 
of the Pade-Borel summation technique. The results obtained iiP^J essentially differ 
from the results evaluated by a double e, 5 - expansion in Ref.'^ The comparison of 
calculated the exponent v values and ratio 2/ a showed the violation of the relation 
z/ = 2/a, supposed irP as exact. 

Ballesteros and ParisP^ have studied by Monte Carlo means the critical behav- 
ior in equilibrium of the 3D site-diluted Ising model with LR spatially correlated 
disorder, in the o = 2 case corresponding to linear defects. They have computed the 
critical exponents of these systems with the use of the finite-size scaling techniques 
and found that a v value is compatible with the analytical predictions u = 2/a. 

In our papei'^ the integrated Monte Carlo simulations of the short-time dy- 
namic behavior have been carried out for 3D Ising and XY models with LR-correlated 
disorder at criticality, in the case corresponding to linear defects. Both static and 
dynamic critical exponents were determined for systems starting separately from 
ordered and disordered initial states. The obtained values of the exponents are in 
good agreement with results of the field-theoretic description of the critical behavior 
of these models in the two-loop approximation.l^^ 

Also, in paper® authors have performed extensive Monte Carlo simulations of 
critical statics of 3D Ising model with LR-correlated disorder with linear defects. The 
Swendsen-Wang algorithm was used alongside with a histogram reweighting tech- 
nique and the finite-size scaling analysis to evaluate the values of critical exponents. 
It was shown that obtained estimates for exponents differ from both previous numeri- 
cal simulation^l^'fi^ and results of renormalization-group (RG) calculations.'^ '^^'^^ 

The present paper is devoted to numerically investigations of critical dynamics 
in short-time regime of 3D Heisenberg magnets with LR-correlated defects. Insertion 
disorder with LR correlations or extended defects must modify a conventional critical 
dynamics of pure Heisenberg ferro- or antiferromagnets which is described by the 
model J and the model G in the classification of Hohenberg and HalperirP^ and 
can lead to relaxational critical dynamics described by the model A with number 
components n = 3 for order parameter. The results of Monte Carlo study of the 
nonequilibrium behavior in the planar magnetics described by 2D XY-model with 
quenched structural defects can be evidence of it.'^D A significant changes in the 
time dependence of the autocorrelation function have been observed in the low- 
temperature phase due to localization of the spin excitations on structural defects. 

The method of short-time critical dynamicJ^ gives the possibility to deter- 
mine both static and dynamic critical exponents modified by LR correlations of 
impurities.'^ In the following section, we introduce the 3D Heisenberg model with 
isotropic distributed linear defects and scaling relations for the short-time critical 
dynamics. In Sec. [3l we derive the critical relaxation in short-time regime for dis- 
ordered Heisenberg systems starting from an ordered initial state and cite values of 
static and dynamic critical exponents obtained with the use of the leading corrections 
to scaling. The final section contains analysis of the main results, their comparison 
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with results of other investigations, and our conclusions. 

§2. Description of the model and methods 

We have considered the following 3D site-diluted ferromagnetic Heisenberg model 
Hamiltonian defined in a cubic lattice of linear size L with periodic boundary con- 
ditions: 

H = -JY^P^pASj, (2-1) 

where Si = {Sf, , 5f ), the sum is extended to the nearest neighbors, J > is the 
short-range exchange interaction between spins Si, and the pi are quenched random 
variables {pi = 1, when the site i is occupied by spin, and pi = 0, when the site is 
empty), with LR spatial correlation. An actual pi set will be called a sample from 
now on. We have studied the next way to introduce the correlation between the pi 
variables for WH model with a = 2, corresponding to linear defects. We start with a 
filled cubic lattice and remove lines of spins until we get the fixed spin concentration 
p in the sample. We remove lines along the coordinate axes only to preserve the 
lattice symmetries and equalize the probability of removal for all the lattice points. 
This model was referred in^''- as the model with non-Gaussian distribution noise and 
characterized by the isotropic impurity-impurity pair correlation function decays for 
large r as g{r) ~ 

In this paper we have investigated disordered Heisenberg magnets with the spin 
concentrations p = 0.8. We have considered the cubic lattices with linear size L = 
128. The Metropolis algorithm has been used in simulations. We consider only 
the dynamic evolution of systems described by the model A in the classification of 
Hohenberg and Halperin.123 The Metropolis Monte Carlo scheme of simulation with 
the dynamics of a single-spin flips reflects the dynamics of model A and enables us to 
compare obtained critical exponents z to the results of RG description of the critical 
dynamics of this model.l^ 

According to the argument of Janssen et alW^ obtained with the RG method 
and e-expansion, one may expect a generalized scaling relation for the A;th moment 
the magnetization 

m(^) (t,T,L,mo) = ^-'^^/^ 

X m^'^) (b'H, b^/^T, b-^L, 6^'«mo) 

is realized after a time scale tmic which is large enough in a microscopic sense but still 
very small in a macroscopic sense. In Eq. ()2-2p . 6 is a spatial rescaling factor, /3 and u 
are the well-known static critical exponents, and z is the dynamic exponent, while the 
new independent exponent xq is the scaling dimension of the initial magnetization 
mo and t = {T — T^jT^ is the reduced temperature. 

The short-time dynamic method in part of critical evolution description of 
system starting from the ordered initial state is essentially the same as the non- 
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equilibrium relaxation method proposed by Ito in Ref.'^''^ for non-equilibrium 
critical behavior study. 

Since the system is in the early stage of the evolution the correlation length 
is still small and finite size problems are nearly absent. Therefore, we generally 
consider L large enough (L = 128) and skip this argument. We measured the time 
evolution of the magnetization determined as follows: 



m{t) 



1/2 



(2-3) 



where angle brackets denote the statistical averaging, the square brackets are for 
averaging over the different impurity configurations, and = pL^ is a number of 
spins in the lattice. 

The question arises how a completely ordered initial state with niQ = 1 evolves, 
when heated up suddenly to the critical temperature. In the scaling form ()2-2p . one 
can skip besides L, also the argument mo = 1, 



m 



(k) 



{t, t) = 6-'=^/^m('=) (^b-H, b^l^T^ . (2-4) 



The system is simulated numerically by starting with a completely ordered state, 
whose evaluation is measured at or near the critical temperature. The quantities 
measured are m{t) and rrS'^\t). With b = t^^^ , one avoids the main t dependence in 
m^^\t) and for k = \ one has 

r) = t-^/^^m(l, t^/'^'T) (2-5) 
= r^l"^ [l + at^l^^T + 0(r2)) . 

For r = 0, the magnetization decays by a power law m{t) ~ t~^/'^^ . If r 7^ 0, the 
power law behavior is modified by the scaling function f^/'^^r). From this fact, 
the critical temperature Tc and the critical exponent f3/vz can be determined. 
The scaling form of magnetization in Eq. ()2-5p is presented as follows: 



In m(t,r) = Int + lnm(l, t^/^'^r) (2-6) 

after differentiation with respect to r gives the power law of time dependence for the 
logarithmic derivative of the magnetization in the following form: 

9.1nm(t,r)|,^o~ti/-, (2-7) 

which allows to determine the ratio 1/uz. On the basis of the magnetization and its 
second moment, the cumulant 

U2{t) = -^-lr^t'^/' (2-8) 

is defined. From its slope, one can directly measure the dynamic exponent z. Con- 
sequently, from an investigation of the system relaxation from ordered initial state 
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Fig. 1. Time evolution of the magnetization m(t) for different values of the temperature T (a) and 
for critical temperature T = Tc = 1.197(2) (b). 

with mo = 1, the dynamic exponent z and the static exponents f3 and can be 
determined. 

The critical dynamic exponent z of the model can be obtained also from the 
time evolution of the raticP^ 



m 



(2) 



{t)\ 



mo=0 



t 



{d-2p/v)/z 



[m{t)r 



-213/vz 



.d/z_ 



(2-9) 



mo = l 



§3. Measurements of the critical exponents for 3D Heisenberg model 

with linear defects 



We have performed simulations on three-dimensional cubic lattices with linear 
size L = 128, starting from an ordered initial state. We would like to mention that 
measurements starting from a completely ordered state with the spins oriented in 
the same direction (mo = 1) are more favorable, since they are much less affected by 
fluctuations, because the quantities measured are rather big in contrast to those from 
a random start with with zero or small initial magnetization (mo ^ 1). Therefore, 
for careful determination of the critical temperature and critical exponents for 3D 
Heisenberg model with linear defects, we investigate the relaxation of this model 
from a completely ordered initial state. 

Initial configurations for systems with the spin concentration p = 0.8 and ran- 
domly distributed quenched linear defects were generated numerically. Starting from 
those initial configurations, the system was updated with Metropolis algorithm. Sim- 
ulation have been performed up to t = 1000 Monte Carlo steps per spin (MCS/s). 

We measured the time evolution of the magnetization m(t) and the second mo- 
ment m^'^\t), which also allow to calculate the time-dependent cumulant U2{t) in 
Eq. (IMD- 

In Fig. [T^ the magnetization m(t) for samples with linear size L = 128 at 
T = 1.191, 1.194, 1.197, 1.20, and 1.203 is plotted in log-log scale. The resulting 
curves in Fig. [1] have been obtained by averaging over 2800 samples with different 
linear defects configurations with 25 runs for each sample. We have determined the 
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Fig. 2. Cumulant U4{T, L) (a) and ratio ^/L (b) as a function of T for lattices with different sizes 
L. 



critical temperature Tc = 1.197(2) from best fitting of these curves by power law. 
The magnetization m{t) at the critical temperature T = Tc is plotted in Fig. [T]3. 

In order to check-up the critical temperature value independently, we have car- 
ried out in equilibrium the calculation of cumulant U4, defined as 



[m 



and the correlation length!^ 

1 



2 sin (tt/L) 



I- 



X 
F 



1 

'Ns 

1 

1 

3 



mi 



E E 

k={x,y,z} n=l 



PjS^ exp (^— ^ 



where {xij,X2,j,xsj) are coordinates of j'-th site of lattice. 
The cumulant Ui{L,T) has a scaling form 



(3-1) 



(3-2) 
(3-3) 
(3-4) 

(3-5) 



(3-6) 



The scaling dependence of the cumulant makes it possible to determine the critical 
temperature Tc from the coordinate of the points of intersections of the curves speci- 
fying the temperature dependence Ui{L, T) for different L. In Fig. [5^ the computed 
curves of U4{L, T) are presented for lattices with sizes L from 32 to 128. As a result 
it was determined that the critical temperature is Tc = 1.197(2). In this case for 
simulations we have used the Wolff single-cluster algorithm with elementary MCS/s 
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Fig. 3. Time evolution of the cumulant Fig. 4. Time evolution of the logarith- 
F2(t) is plotted on a log-log scale at mic derivative of the magnetization 

T — Tc(p). dr him{t, T)|^^y with respect to r is plot- 

ted on a log-log scale. 



step as 5 cluster flips. We discard 256 MCS/s for equilibration and then measure 
after every MCS with averaging over 2048 MCS/s. The results have been averaged 
over 1000 different samples with 25 runs for each sample. 

The crossing of was introduced as a convenient method for calculating of Tc 
inW^ In Fig. [Jb the computed curves of temperature dependence of ratio ^/L are 
presented for lattices with the same sizes, the coordinate of the points of intersections 
of which also gives the critical temperature Tc = 1.198(5). The value of Tc = 1.197(2) 
we selected as the best for subsequent investigations of the Heisenberg model with 
linear defects and with spin concentration p = 0.8. 

In Fig. [3] the cumulant F2{t) and in Fig. [H the logarithmic derivative of the 
magnetization drl^m{t,T)\^^Q with respect to r are plotted on a log-log scale at 
T = Tc- The drlnm{t,T)\^^Q have been obtained from a quadratic interpolation 
between the three curves of time evolution of the magnetization in Fig. [T] for the 
temperatures T = Tc, T = Tc ^ 0.003 and taken at the critical temperature Tc = 
1.197(2). The resulting curves in Figs. [1] and [3] have been obtained at the critical 
temperature by averaging over 3800 samples with 25 runs for each sample. 

In contrast to short-time dynamics of the pure systems we can observe the 
crossover from dynamics of the pure system on early times of the magnetization 
evolution from t ~ 15 up to t ~ 35 MCS/s to dynamics of the disordered system 
with the influence of long-range correlated defects for t > 80 MCS/s. The same 
crossover phenomena were observed in the Monte Carlo simulated critical short- 
time dynamic behavior of diluted 3D Ising systems with point-like defect^^ and in 
behavior of 3D Ising and XY systems with linear defects.'^ 

The existence of different regimes in short-time dynamic evolution of disordered 
system can be demonstrated clearly on critical time behavior of the ratio F2{t) 
(Fig. [3]). We have analyzed the time dependence of the F2{t) in the time interval 
t G [15,35] MCS/s where the F2{t) is best fitted by power law with the exponent 
d/z = 1.464(22). This value of d/z gives the dynamic exponent z = 2.049(31), 
corresponding to the pure 0(n = 3) Heisenberg model.'^D An analysis of the F2{t) 
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Fig. 5. Dependence of the mean square error as a function of the right border of the time interval 

t G [80, tright] 

slope measured in the interval t G [80,300] MCS/s shows that the exponent d/z = 
1.217(3) which gives z = 2.465(6). The dependence of the mean square error cr^ 
as a function of the right border of the time interval t € [80, tright] is presented in 
Fig. [5l The time interval [80, tright] for tright = 300 gives the minimum of errors for 
exponent z in comparison with value z = 2.591(10), which evaluated in the time 
interval t G [80,1000]. 

The same analysis of time dependences of the magnetization m{t) and the log- 
arithmic derivative of the magnetization 9,- lnm(t, t)|^^q in the initial time interval 
gives the exponent values (3/i^z = 0.249(1) and 1/i^z = 0.692(15) with the use 
o{ z = 2.049(31) we obtain the values of exponents u = 0.705(26), f3 = 0.360(9) 
and = 0.510(10). These values are in a good agreement with exponents v = 
0.7048(30), /3 = 0.3636(45) and 13 /v = 0.5158(102) obtained in YleiW^ on the base 
of high resolution Monte Carlo study of pure 3D Heisenberg model critical behavior. 

An analysis of the U2{t) slope measured in the interval t € [80,300] MCS/s 
shows that the exponent d/z = 1.170(136) which gives z = 2.564(368). The value of 
the dynamic exponent z obtained from the ratio F2 in Eq. (j2-9p is more preferred 
because evolution oi F2{t) is less fluctuated then U2{t), and, therefore, a greatly large 
statistics is necessary for obtaining the same quality results from measurements of 
U2{t) as is in the case with F2{t). So, the slope of magnetization m(t) and its 
derivative d\nm{t)\^^Q measured in the interval t G [80,300] MCS/s provides the 
ratio of exponents /3/i/z = 0.150(1) and 1/uz = 0.483(22) which give /3 = 0.311(16) 
and ly = 0.840(47). 

These values of exponents can be compared with results of the field-theoretic 
description of the critical behavior of Heisenberg model with LR-correlated defects 
in the two-loop approximation z = 2.264, v = 0.798, and /3 = 0.384 calculated 
m 

Ref.'i^) for 

case of Heisenberg system with linear defects when the correlation 
parameter a = 2. The some numerical differences of these values exceeding the 
limits of statistical errors of simulation and numerical approximations should not be 
discouraged since the obtained effective values of exponents cannot be considered as 
final. 
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Fig. 6. Dependence of the mean-square errors a of the fits for (a) the cumulant F2, (b) magneti- 
zation m, and (c) the logarithmic derivative of the magnetization dr \nm{t,T) as a function of 
the exponents d/z, P/vz and l/uz for uj/z = (a;/2;)min and time interval [80,300]. 
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Fig. 7. Dependence of global mean-square error (a) (b) ^^/,yz and (c) zii/^^^ for all time 

intervals as a function of the exponent uj/z. 

In the next stage in order to obtain accurate values of the critical exponents, we 
have considered the influence of a leading corrections to the scaling on asymptotic 
values of exponents. We have applied the following expression for the observable 
X{t): 

X{t)=A,t\l + B,t-^/'), (3-7) 

where w is an exponent of the leading corrections to scaling, and are fitting 
parameters, and an exponent 5 = —(3/vz when X = m{t), 5 = d/z when X = 
F2{t) ox X = U2{t), and 6 = Xjvz when X = drlxim{t,T)\^^Q. The expression in 
Eq. ()3-7p reflects the scaling transformation in the critical range of time-dependent 
corrections to scaling in the form of t^^^^ to the usual form of corrections to scaling 
T^'^ in equilibrium state for time t comparable to the order parameter relaxation 
time tr ~ C'^ikC)^ 

We have used the least-squares method for the best approximation of the simu- 
lation data X{t) by the expression in Eq. (j3-7p . Minimum of the mean square errors 
as of this fltting procedure determines the exponents 5 and u/z. As example, we 
plot in Fig. [6] for time interval [80,300] the a^/^ (a) for the cumulant F2 as a func- 
tion of the exponent d/z for u/z = 0.40, the u^/jy^ (b) for the magnetization m as a 
function of the exponent /^/I'z for u/z = 0.24 and the 0"i/^^ (c) for the logarithmic 
derivative of the magnetization as a function of the exponent l/i^z for uj/z = 0.41. 
All obtained minimal values of 6 are averaged for all time intervals. The distribution 
function p{5) for different time intervals is presented in Fig. [8] for instance 5 = d/z. 
The global mean-square error 

^5 = VE^Ii (S,-5r-Au/j:^^,AU (3-8) 
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Fig. 8. The distribution function p{d/z) for different time intervals. Solid line is corresponds to 
calculated averaged value d/z — 1.329 
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Fig. 9. Dependence of global mean-square error A/jfjy^ as a function of the exponent uj/z for dif- 
ferent values of T. 

is calculated on base of averaged values 5, where N is number of time intervals, which 
demonstrates minimum of errors. The dependence of global mean-square error A^/^ 
(a) for the cumulant F2, A^j^^ (b) for the magnetization m and Aij^^^ (c) for the 
logarithmic derivative of the magnetization as a functions of the exponent uj / z are 
presented in Fig. [71 The critical exponents are calculated for w/z which corresponds 
to minimum of A^. 
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Table I. Values of the obtained critical exponents and comparison with results of renormalization 
group (RG) and Monte Carlo (MC) calculations 





Z P/v V P LO 


mo — 1, LR system stage 
Too = 1, pure system stage 


2.257(61) 0.510(78) 0.770(74) 0.393(77) 0.786(45) 
2.049(31) 0.510(10) 0.705(26) 0.360(9) 


LR system (a — 2) 


Prudnikov, et al, 2000, (Ref.^^') RG d = 3 
Prudnikov, et al, 2010, (Ref.^"*') RG d = 3 
Blavats'ka, et al, 2001, (Ref.^"') RG d = 3 


2.264 0.482 0.798 0.384 
2.291(29) 0.490(5) 0.766(17) 0.375(5) 

0.88 


Pure system 


Prudnikov, et al, 2008, (Ref!^) RG e-exp. 
Guida, et al, 1998, (RefP) RG d = 3 

RG e-exp. 
Chen, et al, 1993, (RefP) MC 


2.020(7) 

0.5178(13) 0.7073(35) 0.3662(27) 
0.5188(23) 0.7045(55) 0.3655(45) 
0.5158(102) 0.7048(30) 0.3636(45) 



The value of exponent d/ z = 1.336(35) was computed with oo/ z = 0.40(1), value 
fi/vz = 0.228(7) was computed with u/z = 0.24(1) and value l/uz = 0.589(4) 
was computed with uj/ z = 0.41(1). It can be determined the values of the critical 
exponents z = 2.245(60), v = 0.757(26), 13 = 0.388(15). For the averaged value 
uj/z = 0.36(7) it was computed values of the exponents d/z = 1.329(36), j3/vz = 
0.226(29), and 1/vz = 0.575(40). On the base of these values, we determine the 
final values of the critical exponents z = 2.257(61), v = 0.770(74), /3 = 0.393(77), 
and w = 0.786(45). The statistical errors for exponents are estimated by dividing all 
data into five sets. 

The dependences of global mean-square error Apj^z as a function of the exponent 
UJ / z for different values of T are presented in Fig. and they demonstrate that 
temperature T = 1.197, which was choosen as critical, gives the minimal value of 
fitting errors. 

§4. Analysis of results and conclusions 

The present results of Monte Carlo investigations allow us to recognize that the 
short-time dynamics method is reliable for the study of the critical behavior of the 
systems with LR-correlated disorder and is the alternative to traditional Monte Carlo 
methods. But in contrast to studies of the critical behavior of the pure systems by the 
short-time dynamics method, in case of the systems with quenched disorder corre- 
sponding to randomly distributed linear defects after the microscopic time tmic — 10 
there exist three stages of dynamic evolution. In the time interval of 15 — 35 MCS/s, 
the power-law dependences are observed in the critical point for the magnetization 
m{t), the logarithmic derivative of the magnetization drl'nrn{t,T)\^^Q, and the cu- 
mulant F2{t), which are similar to that in the pure system. In the time interval 
[80,300], the power-law dependences are observed in the critical point which are de- 
termined by the influence of disorder. In the intermediate time interval the crossover 
behavior is observed in the dynamic evolution of the system. However, careful anal- 
ysis of the slopes for dynamical characteristics reveals that a correction to scaling 
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Fig. 10. The cumulant f/4 as a function Fig. 11. Dependence /(C/4) = [in C/4/ ln(T - Tc)]"^ 
of a variable {T -T,)L^^'' {l + bL--^) as a function of 1/L. 

( • is corresponds to T = 1.190, 
★ - r = 1.195, T - r = 1.200, 
■ - T = 1.205, A - T = 1.210). 



should be considered in order to obtain accurate results. 

The cumulant U4 as a function of the scaling variable (T—Tc)L^^'^ {l+bL~'^) with 
using correction to scaling procedure for calculated values u = 0.770 and uj = 0.786 
is shown in Fig. [TUl Approximation of function /([/4,1/L) = [lnC74/ln(r — Tc)]^^ 
(Fig. [TT]l gives possibility to estimate the exponent = 0.758(10), which is in a good 
agreement with final value u = 0.770(74) (Table |I|). 

The dynamic and static critical exponents were computed with the use of the 
leading corrections to scaling for the 3D Heisenberg model with linear defects and 
their values z = 2.257(61), u = 0.770(74), (3 = 0.393(77), and w = 0.786(45) can 
be considered as final. In a summary Table lU we present the values of critical 
exponents z, ji/v, v, /3, and oj obtained in this paper by comprehensive Monte Carlo 
simulations of the short-time critical evolution of the diluted 3D Heisenberg model 
with linear defects from an ordered initial state with mo = 1. For comparison, we 
give in Table |T] the results of renormalization group (RG) and Monte Carlo (MC) 
calculations of these exponents for pure 3D Heisenberg modet^"^''^^ and diluted 
3D Heisenberg model with linear defects .'^''^''^^ 

The obtained in this article values of exponents demonstrate very good agree- 
ment in the limits of statistical errors of simulation and numerical approximations 
with results of the RG field-theoretic description from Ref.,'^ calculated with the 
use of the Pade-Borel (PB) summation technique to d = 3 expansion series, and 
particularly with results from Ref.,!^ where the Pade-Borel-Leroy (PBL) and the 
self-similar approximation (SSA)PSli resummation methods were also applied to series 
from Ref."^^-' The obtained value of the correction-to-scaling exponent w = 0.78(31) 
demonstrates a sufficiently good agreement with value of uo = 0.88, obtained in 
Ref.'SSJ by the RG field-theoretical method with fixed dimension {d = 3) for Heisen- 
berg model with isotropic distributed linear defects (a = 2). 

The obtained results confirm the strong infiuence of LR-correlated quenched 
defects on the critical behavior of the systems described by the many-component 
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order parameter. As a result, a wider class of disordered systems, not only the 
three-dimensional diluted Ising model, can be characterized by a new type of criti- 
cal behavior induced by randomly distributed quenched defects and effects of their 
spatial correlations. 
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